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T^ij- , Abstract 

'^ \ We study deeply virtual Compton scattering on a virtual pion that is emitted by a proton. 

C^ ' Using a range of models for the generalized parton distributions of the pion, we evaluate the 

cross section, as well as the beam spin and beam charge asymmetries in the leading-twist 

f^ i approximation. Studying Compton scattering on the pion in suitable kinematics puts high 

demands on both beam energy and luminosity, and we find that the corresponding require- 
ments will first be met after the energy upgrade at Jefferson Laboratory. As a by-product 
of our study, we construct a parameterization of pion generalized parton distributions that 
C^ I has a non-trivial interplay between the x and t dependence and is in good agreement with 

form factor data and lattice calculations. 



1 Introduction 

The concept of generalized parton distributions (GPDs) is a versatile tool to describe hadron structure 
at the quark-gluon level and has given rise to vigorous theoretical and experimental activities. Among 
the attractive features of GPDs are the possibilities to connect ordinary parton densities with elastic 
form factors [1] and to explore the spatial distributions of partons inside a hadron [2,3]. Reviews of 
this extensive subject can be found in [4-7]. 

The pion plays a special role in the low-energy sector of QCD as the lightest bound state and the 
Pseudo-Goldstone boson associated with chiral symmetry breaking. Given the difficulty to perform 
high-energy experiments with a pion in the initial state, our knowledge of its internal structure is, 
however, scarce compared with what is known about the nucleon. Currently, the principal sources of 
information are the spacelike electromagnetic form factor of the pion [8-11] and its parton densities 
extracted from Drell-Yan production with pion beams [12-14]. Measurements constraining the GPDs 
of the pion would provide a natural extension of this knowledge. On the theoretical side, the pion 
GPDs have been studied in a number of dynamical models [15] and on the lattice [16-19]. Important 
theoretical investigations have been performed for a pion target in the first instance, because of its 
relative simplicity as a spin-zero particle, see for instance [20,21]. 

The purpose of the present work is to estimate how pion GPDs may be investigated in deeply 
virtual Compton scattering (DVCS), which is the theoretically cleanest and most advanced among 
the hard exclusive processes involving generalized parton distributions. We consider the reaction 
ep — > e'jTT'^n at small invariant momentum transfer between the proton and neutron. In the one-pion 
exchange approximation, the reaction is then described by the emission of a slightly off-shell pion from 
the proton, followed by the scattering process evr"'' — > e-fn^. This can be seen as an analog of the 
reaction ep — > e7r"^n, which has been used to extract the electromagnetic pion form factor for all but 
the smallest values of momentum transfer [9-11]. Two mechanisms contribute to evr — > cjn, namely 
virtual Compton scattering and the Bethe-Heitler process, as shown in Fig. [TJ Suitable strategies for 
isolating the Compton signal are the same as those for scattering on a nucleon target, which have been 
elaborated in detail [22] and successfully used in experiment [23,24]. The corresponding expressions 
for the pion can be found in [25], where also numerical estimates for evr — > ejn have been given. 
We note that in suitable kinematics, the reaction ep -^ e^TT~^n may also be used to study virtual 
Compton scattering on the pion in the backward region, whose description involves the so-called 
transition distribution amplitude from a pion to a photon [26]. In experiments with a real photon 
beam, the lepton-pair production process 7p — > e+e~7r"*"n can provide access to timelike Compton 
scattering 77r — > 7*7r on the pion, which is closely related to DVCS by crossing [27]. 

Our paper is organized as follows. In the next section we present in some detail the kinematics of 
ep -^ e'-fir^n. In Sect. [3] we give the basic equations for the analysis of this process in the one-pion 
exchange picture and in the framework of generalized parton distributions. We also briefly discuss 
the validity of the one-pion exchange approximation. In Sect. [5] we present a model for the GPDs of 
the pion, paying special attention to their t-dependence. In Sect. Owe finally give estimates for cross 
sections and asymmetries in the HERMES experiment and at Jefferson Lab, before summarizing our 
main findings in Sect. [H 

2 Kinematics 

In this section we discuss the kinematics of the reaction 

e(/) + p(p) ^ e(/') + 7(^0 + vr+(p;) + n(p') , (1) 





Figure 1: Graphs for ep — > e^7r~^n in the one-pion exchange approximation. Contributing subprocesses 
are virtual Conipton scattering on a pion (left) and the Bethe-Heitler process (right). The crossed 
Bethe-Heitler graph (not shown) has the photons attached to the lepton line in opposite order. The 
blob marked with F.,^ represents the electromagnetic pion form factor. 



with four-momenta given in parentheses. We write 

q = l-l', Pn=p-p' (2) 

for the four-momenta of the virtual photon and the virtual pion, and use the standard variables 

Q' = -q\ W' = ip + qf, s = {p + l)\ XB = ^, y = ^ (3) 



2p ■ q p ■ I 



for deep inelastic scattering processes. The variables 



t = {p-p')\ x^ = ^ (4) 

p ■ L 

describe the emission of the virtual pion from the proton target, where x,r is the fraction of energy 
that the virtual pion takes away from the proton in the ep cm. The azimuthal angles of the final-state 
electron and neutron in that frame are respectively denoted by ij^e and ipn, where the z-axis is chosen 
along the lepton beam momentum. We write ?7i/v and 171-,^ for the nucleon and pion masses, and 
neglect the lepton mass throughout our work. An important role is played by the kinematic limit 

-i>-to = f^, (5) 

which is readily obtained from the expression p'^ = — t(l — x^) — x^rn^ for the squared transverse 
momentum of the neutron in the ep cm. We need two more variables to describe the 7 and vr in the 
final state, namely 

K = {p.-p'.? (6) 

and the azimuthal angle (p-,^ between the plane spanned by p^ and p'^ and the plane spanned by I and 
I' in the cm. of vrx in the final state. For the sign of (pT^ we follow the usual convention for DVCS 
on a proton target |1| For our later discussion it is useful to introduce further variables, which refer to 
the subprocess evr — > e^ir on the virtual pion target, namely 

Sn = [Pn + q) , Xb = , y-K = 7 • (7) 

2p-K-q Ptt -i 



One thus obtains (pn from the angle ph in [28] by replacing P ^ pn and Ph 



One finds 

Stt = Q^(^^ - 1) - 2cos(V'e -tpn) \/(l - y)Q'^ - {yxsmNY \/(l - x^)(to - 1) 
\Xb / 

+ 2yxBXT,m% + (1 - yxB)t . (8) 

To select kinematics where DVCS on a pion can be described in terms of generalized parton 
distributions, we take the Bjorken limit 

Q ^ oo at fixed y, xb, Xt^, t, t,,-. (9) 

The squared cm. energies s, W"^, and s-,^ then become large, whereas the squared momentum transfers 
t and i,v are small compared with Q^. The relation ([8]) then implies 

X77- 





y-K 


= ?/ 

X-TT^B 


x^2/ R 


^ s. + Q2 





and 

(11) 

with corrections of order m/Q, where m? represents the small scales m'^, m^, t and iyri whose 
magnitudes we do not distinguish at this point. The phase space element of the process ([T]) can be 
written as 

The interpretation of the process ([1]) in terms of DVCS on a virtual pion target puts several 
conditions on the kinematics, which we now discuss. First of all we impose an upper cutoff on |t|, 

I''! S I'^lmax ) l^'JJ 

to ensure that the p — > n transition is dominated by virtual pion emission. Since |t|max rnust be bigger 
than —to in ([5]), this implies a maximum value for x^, 



Ir 



Xtj- ^ 3^7rmax ^ 



V r^ + 4r — r 



with T = %^ . (14) 



m^ 



Since we want the subprocess 7*7r — > jn to be in Bjorken kinematics, we further impose lower cutoffs 
According to (jlip this implies lower limits on x^r and on y, 

, ymin ^ , \^^) 



where ymax is an upper limit on y we will later impose both for theoretical and for experimental 
reasons (see Sect. [5]). The relation ([TT]) also restricts the largest possible values of Q^ to 

^max ^^ ■'■^Trmax "max ''Trmin ' \^ ' ) 

For a clean physical interpretation of the reaction ([T]) as DVCS on a weakly off-shell pion target, it is 
desirable to take rather small |t|max and rather large s^j^^jj^ and Q^j^. With (fH|) and the first relation 
in (I16p . this only leaves enough phase space for x-j^ if the total cm. energy s is sufficiently large. 



3 Calculation of the cross section 

In this section we give some basic formulae for the process ep — > e^ynn in the one-pion exchange 
approximation and discuss the vahdity of this approximation. 

3.1 The one-pion exchange approximation 

In the one-pion exchange approximation, the amphtudes for ep — > e^nn and for evr — > e^n are related 
as 



V^girNN r-u\ A^ T-r^\ ^^ 



m 



2 



Mep = U{p'h5u{p) "^r_^ ^(*) -^- ' ^(*) = ^2_/ (18) 

with the pion-nucleon coupling QnNN = 13.05 [29]. Here we have introduced a phenomenological 
factor F{t) to soften the pion-nucleon vertex when the pion virtuality t becomes large compared to 
m^. In our calculations we will take A = 800 MeV from [30]. The ep cross section is then given by 

d^a{ep -^ e^irn) 1 



dy dQ^ dipe dU dcj)^ dt dx^ d^„ 128 (27r)8 (s - m^) y/{s^ + Q'^ + ty - is^t 

spins 

where X^spins sums over the polarizations of the final-state electron and photon. In (I19p we have 
further averaged over the polarization of the proton target and summed over the polarization of the 
outgoing neutron, but kept the lepton beam polarization fixed. Defining the cross section 

dy^dQ^dt^dcP^ 32(27r)4x^(s-m2,)V(s^ + Q2 + t)2_4s^t ^^J ^"' 

on a virtual pion target, we have the simple relation 

d'^a{ep^e-fTTn) _ qInn ij:^u\l^ ~^ d'^ajeTT ^ e-fir) 



dy dQ"^ dtj^ d(l)TT dt dxj^ ^ Svr^ {rn^ — tY dy.,^ dQ"^ dtj, 

where we have integrated over the angles ipe and f/'n- 

For a rough estimate one may neglect the dependence of the evr cross section (j20p on the pion 
virtuality t. Integrating the factor of proportionality in ()2ip over i, one then obtains 

d'^a{ep ^ e-iTTn) ^ f _(iV(e7i^^^^77r)_ 



dy dQ"^ dt-jt d(j)T^ J ' dyT^ dQ"^ dtj^ d^Tr 

with 

n{x^,\tU.) = xJ-^ I dt[F{t)]'^-^^, (23) 

1 1 1 max 

where to(a^7r) is given in ([5]). Note that the evr cross section on the r.h.s. of (f22]) depends on Xt^ via 
s,r ~ x^^y^^s — Q^. With the form factor of the pion-nucleon vertex taken in (fT8]l. the integral in (p3]) 
can be done analytically using 

^^[^(*)]'r^^ = ^^i-4^-#7-^- (24) 

'- -I {ml - t)2 A2 - m2 ml-t A^ - t ml - t 
In Fig. [2] we show the pion fiux factor Il{xTr, |i|max) as a function of Xjr for two values of |t|max- 







Figure 2: The pion flux factor Il{x.j^, |t|max) defined in (f23]) . shown for |t|max = 0.3 GeV^ (solid) and 
|i|max = 0.5 GeV^ (dashed). 



3.2 Validity of the one-pion exchange approximation 

The validity of the one-pion exchange approximation in the process ([T]) cannot be taken for granted, 
and it is natural to see what is known for similar processes. 

As already mentioned in the introduction, an important process in this context is ep -^ evrn, where 
the one-pion exchange approximation yields the subprocess eir — > evr, from which the electromagnetic 
pion form factor can be extracted. The most recent measurements [10, 11], as well as previous data, 
provide a clear indication that pion exchange cannot be the only contribution to this reaction. This 
is because the cross section ctt for ^*p — > vrn with transverse 7* polarization in the ■~f*p cm. is seen to 
be clearly nonzero, although at low |f| it is smaller than the cross section cjl for a longitudinal 7*. To 
understand the implications of this observation, let us take the case |t| = |t|min) where by definition 
the angle between the proton and neutron momenta is zero in the j*p cm. Due to angular momentum 
conservation, the subprocess 7*7r -^ vr can then only proceed for longitudinal photon polarization, i.e., 
pion exchange contributes only to 0"/,. In terms of t-channel exchanges, the presence of (Tt is thus an 
indicator for the exchange of states with nonzero spin, with the p being an obvious candidate because 
of its relatively low mass. A corresponding model calculation in [31] indeed yields a nonvanishing ut, 
although it undershoots the measured values of [10] and [11], which were taken at j*p cm. energies 
of W = 1.95 GeV and W = 2.22 GeV, respectively. In [10] this mismatch was ascribed to possible 
resonance contributions in the vrn channel. At |t| = |i|min there is no contribution to ul from p 
exchange because of parity conservation in the subprocess j*p — > vr, so that for the extraction of the 
pion form factor a separation of ai and ax should considerably enhance the contribution from pure 
pion exchange. 

We note that the situation in 7*p -^ ^im is different. Even at |t| = |i|min both subprocesses 
7*7r — > 77r and 7*p — > 77r can take place for transverse 7* polarization, which is dominant in Bjorken 
kinematics according to the factorization theorem for DVCS. Whereas taking low |t| will enhance pion 
exchange also in this case, it may be useful to assess the quantitative importance of 7*/? — > 77r. This 
would require information about the GPDs for the p — > vr transition. Simple helicity counting [32] 
shows that at twist-two level there are two GPDs involving the axial current and one GPD involving 
the vector current for the quarks. Their x moments are accessible to an evaluation in lattice QCD, 
which may thus help to estimate the role of p exchange in the reaction ep — > ejirn. In particular, the 
lowest moment of the vector current GPD gives the electromagnetic p — > vr transition form factor. 
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which also enters the Bethe-Heitler process in ep -^ cyn. A further constraint on the p ^ n transition 
GPDs can be obtained by taking the hmit of soft pion momentum, in analogy to what has been done 
for the A^ — > vr transition in [33]. 

Concerning resonance contributions, we recall that the ^*p cm. energy W is large in Bjorken 
kinematics. In experiments suited to investigate DVCS, W will thus be much larger than in the pion 
form factor measurements [10, 11]. More problematic are possible resonances in the vrn channel of the 
final state, which are of course not taken into account by the one-pion exchange approximation. The 
relevant invariant mass is 

M2„ = (p; + p'f =ml + 2m^E^ +U-t 
1 



m^ + ml-\ 



{tnO -U)il -X,r) -t(l -Xg) + ^ml 

1-% 



+ 2cos(,/.^ + iPe- ^n)yj{l - x^){l - x|)(to " t){Uo - U) +o{^y (25) 

Here Et^ is the energy of the outgoing pion in the proton rest frame, and tvrO is the upper kinematic 
limit of tjj. In the Bjorken limit one has 

Uo^x-^t-^ (26) 

i x^ 

with corrections of order (x^tji)'^ /Q'^. In our numerical calculations we use the exact expression of 
t^jQ, so that any (j).,^ dependent term exactly vanishes at the kinematical point t,r = ^ttOi where (j)-,^ is not 
defined. We see from (|25|) that to avoid resonances contributions in the vrn channel it is advantageous 
to have low x^ and relatively large |t7r| (while still respecting the condition It^^l <C Q^ for Bjorken 
kinematics). 

A different process relevant in our context is deep inelastic scattering with a leading neutron in 
the target hemisphere, ep — > en + X. In the one-pion exchange approximation, this gives access to 
inclusive deep inelastic scattering 7*7r — > X and thus provides information of the parton densities of 
the pion. There is a number of theoretical investigations focusing on very high energies, as achieved 
in the HERA collider experiments [34]. In particular, the studies [35-37] have investigated the role of 
p and also of 02 exchange in the framework of Regge theory. Furthermore, rescattering of the neutron 
has been studied in [37-40] and is typically found to be non-negligible even for Q^ of several GeV . 
Given the high-energy limit underlying these investigations, we find it difficult to assess the situation 
for deeply virtual Compton scattering at significantly lower energies. 

In summary, we find that existing theoretical investigations of similar processes cannot readily be 
used to quantify effects beyond the one-pion exchange approximation for ep -^ e^Trn. They emphasize, 
however, the importance of taking |t| as small as possible. Working at low x^^ will in addition help 
to avoid resonance effects between the outgoing neutron and pion. The incorporation of p exchange 
into the theoretical analysis should be practically feasible (at least at the level of estimates) if one 
could gain some information on the size of the /) — > vr transition GPDs, for instance from lattice 
calculations. 

3.3 Compton scattering on the pion 

Let us now recall the basics of the process evr — > e'jn in the Bjorken limit, which have been elaborated 
in detail in earlier work [25]. The analysis of this reaction proceeds in close analogy to the well-known 
case of a proton target [22]. It is in fact simpler because the pion has spin zero and thus involves 



fewer GPDs and form factors than the proton. Throughout this section, we retain only the leading 
terms in the l/Q expansion, unless explicitly indicated. 

We decompose the amplitude for evr — > e^fii into contributions from Compton scattering and from 
the Bethe-Heitler process. 

Men = Mycs + Mbh ■ (27) 

The corresponding decomposition of the differential cross section for evr — > ejn and thus also for 
ep -^ e'yirn reads 

da = daycs + daBH + dawT , (28) 

where duiNT is the interference term between the Bethe-Heitler and Compton processes. In the 
Bjorken limit we have for the squared Compton amplitude 



^ |2 4e6 1-^^ + ^2/2 
2_^\Mvcs\ =7)2" ~ l^'^l 



spin 



Q' 



yl 



with 



nAC,t^) = Y,4J dxH^{x,ttn) 



^ — X — ie ^ + X — ie 



(29) 



(30) 



at leading order in a^- Here Ht^ is the GPD for quark flavor q in a vr^ as defined in [5], e^ is the 
quark charge in units of the positron charge e, and X^spin sums over the polarization of the final-state 
photon. The skewness variable is given by 



'-s 



The squared Bethe-Heitler amplitude can be written as 

^ 2 We^ l-x^t^o-t^l-y„ + yl/2 

spin 

where the factor 



[FniU] 



P 



-s'u' 



{As — B cos (J)tj)[Au — B cos (/)7r) 



(31) 



(32) 



(33) 



with s' = {I' + q')^ and u' = {I — q')"^ comes from the lepton propagators. This factor is unity in 
the Bjorken limit, but it can deviate quite significantly in experimentally relevant kinematics. Up to 
relative corrections of order x^m? / Q"^ , one has 



As 



Q' 



^1 - y-K)t-K 



y-n 



-^u 



[I - yn)Q'^ - U 
y-K 



B 



2Q 

y-K 



{I - y^){l - xl)(t^o - tn) , (34) 



so that for 4(1 — xj^) (tirO — ^tt) '^ (1 — y-K)Q'^ one can have P close to zero. In our numerical calculations 
we use the exact expressions of s' and u' in (|33p . The interference term between Compton scattering 
and the Bethe-Heitler process reads 

-^VCS-^BH + C-C. - ei —^ — -p^^ -^7r(l7r) 



P 



spin 



^B 



y^^Lfi 



cos (pTT H-e ri-n + Pi ■ sm (p-„: Im n. 



y7r\/l - yiT 



VI -y-n 



(35) 



where e£ = ±1 is the charge of the lepton beam and P^ = ±1 its hehcity. 

For lack of better knowledge, we will ignore the off-shellness of the incoming pion when evaluating 
i^^(t^) and H^{x,^,tTr). In kinematical factors we take, however, the virtuality t of the initial pion 
instead of m^. As can be seen in (j26p . this has an important effect on ^,^0 and thus on the factors 
t-iro — tn in (f32l) and (f35l) . The approximation (i22]) . where the t-dependence of da {en — > e^n) is 
neglected, should therefore be used with caution, especially for small \tT^\. 

In the Bjorken limit, the Bethe-Heitler process dominates over Compton scattering (unless yjr is 
very small). This is because |A1vcsP/|-^BHp ~ ^tt/Q^ according to (p9]) and ([32]) . In this situation, 
privileged access to the Compton amplitude is provided by the interference term (|35]) . which can be 
separated from the cross section by reversing the beam charge ei or the beam helicity Pg. We remark 
that there are also Pg dependent terms in daycs- As can be seen in [25] they are, however, suppressed 
by 1/Q compared with the dominant term given in ()29p . 

To conclude this section we remark on the process en — > e^7r~p, which is accessible through 
incoherent scattering on nuclear targets. Comparing the subprocesses evr^ -^ e'yn^ and evr" — > e77r~, 
we find that the amplitudes for the Bethe-Heitler process have opposite sign, whereas those for the 
Compton process are equal. This is because the 'ymr three-point function changes sign under charge 
conjugation while the ^tttt^ four-point function remains the same, and therefore holds even beyond 
the leading approximation in 1/Q. As a consequence, the relations 



daycsien - 


-^ e77r~p) = 


daycsiep- 


-^ e^ir~^n) , 


daBHien - 


-^ ejir^p) = 


daBuiep - 


> ejTT'^n) , 


d(TiNT(en - 


-> e^7r~p) = 


-do-INT (ep - 


-^ ejTT~^n) 



(36) 

hold as long as the one-pion exchange approximation is valid, whereas they will be invalid if for 
instance interference between tt and p exchange is important. In typical fixed-target kinematics da^n 
is much larger than duiNTi so that the relations ([36]) can principle be tested experimentally: when 
going from ep — > e'jn^n to en — > e^'jr~p, the overall cross section should approximately remain 
the same, whereas the beam spin or beam charge asymmetry should change sign. We note that a 
corresponding consistency check for the one-pion exchange approximation in ep — > e7r~^n and en -^ 
e'K~p was performed in the extraction [11] of the pion form factor. 

4 The generalized quark distribution of the pion 

In this section we describe the model for the pion GPDs which we will use in our numerical studies. 
We generate a dependence on the skewness ^ by the model of Musatov and Radyushkin [41], 



/■I /-i-l/SI 
H^{x,U^)= dp da6{x-p-^a)h{P,a)H^{p,0,U) (37) 



with 



,,,,._ r(25 + 2) [{1 - m)' - a^r (^.. 

MA«)-22mr2(6+i) (i_|/3|)2W ■ (38) 

For the profile parameter we will take either 6 = 2 or 6 = 1. The forward limit of the GPDs is given 
by the parton densities in the pion as H${x,0,0) = q^ix) for x > and H${x,0,0) = —qT^{—x) for 
x < 0. As an input we will take the parameterizations of SMRS [12] or of GRS [13], both at scale 
^ = 2GeV. 




U [GeV" 



Figure 3: Data for the electromagnetic pion form factor from [8-11], compared to a fit specified by 
(jiT]) . (gSD, and the parameters in (03]). 



The simplest way to model the tjr dependence is a factorized ansatz 
which automatically satisfies the sum rule 



^e, / dxH^{x,0,t^) = F^{t^). 



(39) 



(40) 



Both theoretical considerations [42] and lattice QCD calculations [16, 17] indicate, however, that the 
dependence of the GPDs on i^ and x is correlated. As a model ansatz we will take an exponential 
iyr dependence with an x dependent slope. Such an ansatz has proven to be quite successful for the 
proton [43,44]. Following [43] we set 



H^^{x,0,U) = //^(x,0,0) exp[t^/(|x|)] 



with 



1 



fix) = a' {I -xyin- + B{l- xy + Ax{l - xy , 



(41) 



(42) 



where a' = 0.9 GeV~^ is taken in accordance with Regge phenomenology. We fit A and B to describe 
the pion form factor through the sum rule (j40p . using the data [8-11] as selected in [45]. With the 
fitted values 



^ = 2.19GeV-^ 
A = 1.35GeV-2, 



S = -0.38GeV-2 
B= 0.58GeV-2 



for GRS, 
for SMRS 



(43) 



we get a good description for FT^it-,^) with both parameterizations of the parton densities, as shown in 
Fig. [3l An intuitive physical interpretation of the function f{x) is obtained in the impact parameter 
representation [2]: for x > the average of the squared transverse distance between the quark and 
the center of momentum of the pion is {b'^)x = 4/(x). The results of our two fits are rather similar 
and yield physically reasonable values, as shown in Fig. 21 
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Figure 4: The average transverse distance \/{b'^)x between a n-quark and the center of momentum of 
the vr"'', as obtained in the two fits of Fig. [3l 



As a caveat we note that the sum rule (I40p only constrains the valence quark distributions, given 
by ii^^(a;,0, ijr) + H'^{—x,0,tTr) = [qnix) — qTrix)] exp[i^/(x)] for x > 0, and is insensitive to the sea 
quarks. Since sea quarks mix with gluons under evolution, one may expect that the tjr dependence 
in this sector is different from the one for valence distributions. The ansatz (14 ip with a common 
tvr dependence for valence and sea quarks (including the strange sea) may hence be regarded as 
oversimplified. Given, however, that even in the forward limit the sea quark distributions in the pion 
are poorly known at present, we deem this ansatz to be acceptable for our purpose. 

Using (|4T]) and (l42|) we can also evaluate the second moment 



^2o(*7r) 



dxxH^{x,0,tTr) 



(44) 



which has been evaluated in lattice QCD [16,17]. In Fig. [5] we compare the results of our fits (|43p with 
a monopole parameterization of the lattice data given in [17]. Although not perfect, the agreement 
is quite good, and certainly much better than the result of the factorized ansatz (I39p . We note that, 
in our parameterization, the contribution of sea quarks to the moment ()44p is below 30% at t^r = 
and smaller at nonzero t^^. 

Notwithstanding the success of the ansatz given by ()4ip and ()42p in reproducing the lowest two 
moments of the pion GPD, some cautionary remarks from the theoretical side are in order. As 
discussed in [43], the asymptotic behavior of the pion form factor at large negative t in our ansatz 
is controlled by the Feynman mechanism and given by the Drell-Yan relation FT^{t) ~ |t|~'''^+^''^. 
Here /3 describes the behavior q{x) ~ (1 — x)^ of the parton densities in the limit x — >■ 1 and is 
predicted to be /? = 2 for the pion [46,47]. Asymptotically the Feynman mechanism hence gives 
a contribution FT^(t) ~ |t|~'^/^ to the form factor, which is power suppressed compared with the 
contribution FT^{t) ~ |i|~^ from the hard-scattering mechanism [46,48], where it is understood that 
both power laws are modified by logarithms. The ansatz (|4ip . (|42p does not have the form generated 
by the hard-scattering mechanism for H!^{x,£,,t) at large t [49]. However, it turns out that the parton 
density parameterizations [12,13] at /i = 2 GeV both have a large-x power /J ~ 1 within less than 10%, 
in stark contrast with the prediction /3 = 2 from power counting. Thus, our ansatz gives FT^{t) ~ |t|"^ 
at large negative t, which is compatible with the monopole behavior that describes very well the 
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Figure 5: The second moment (H^ of the pion GPD H^, obtained with the two fits of Fig. El 
Also shown is the result of the factorized ansatz ([39]) with the GRS parton densities. The band 
corresponds to a parameterization of lattice results given in [17], ^20(^1-) ~ ^2oi^)/ {^ ~ ^-k/M'^) with 
^fo(O) = 0.261(5) and M = 1.37(7) GeV. 



available data [45]. We shall not pursue this issue further, and use our ansatz for H^{x,^,t-,^) as a 
simple candidate form that is not in contradiction with phenomenological constraints. 



5 Cross section estimates 

In this section, we calculate the cross section for ep -^ en^ir and its dependence on the beam charge 
and beam helicity. Since the unpolarized cross section is dominated by the Bethe-Heitler process, the 
corresponding results are largely model independent (as long as the one-pion exchange approximation 
is adequate). On the other hand, we will see that the charge and polarization asymmetries for the 
lepton beam are quite sensitive to the model we assume for the pion GPD. We will compare four 
simple model scenarios: 

1. a skewness dependence generated by ([37|) and ([38]) with 6 = 2 and a t^ dependence given by 
[1]) to (j43p . with the GRS parton densities in the pion. 



2. the same as model 1, but with the SMRS parton densities, 

3. the same as model 1, but with 6=1 instead of 6 = 2, 

4. the same as model 1 , but with a factorizing t^ dependence ([39]) . We use this model for the sake 
of contrast, although it is disfavored by theoretical considerations and lattice data. 

For the pion form factor we use a monopole parameterization F^(f^) = l/(l — tir/M^) with M = 
714 MeV. This provides a very good description of the experimental data, as shown in [45]. 

As explained in Sect. [21 we impose minimal values for Q^ and s-n- and at the same time a maximal 
value for |t|, which requires a sufficiently large energy. We therefore concentrate on typical kinematics 
for HERMES and for the planned Jefferson Lab upgrade to 11 GeV beam energy. The leading- twist 
interpretation of DVCS demands that |i^| <C Q^, so that we also put a cut |t7r| < I^Trlmax- We 
finally impose a maximum value on y (and thus on y^r)- On the experimental side, this ensures 
that the scattered lepton has sufficient energy to be detected and identified. On the theoretical side. 
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this improves the l/Q expansion underlying the approximate formulae p2l) and ()35p . since there are 
subleading terms in l/Q that come with a factor 1/-^/! — y-,^ relative to the leading terms. An example 
of such a subleading term is found in the propagator factor P, see (I33p and (1340 . 

Let us first consider the case of HERMES, with a beam energy E^, = 27.6 GeV in the proton 
rest frame. We impose a lower limit ^^ > 2.57° on the angle between the momenta of the final- 
state photon and the lepton beam in the target rest frame. This value corresponds to the maximal 
geometric acceptance of the electromagnetic calorimeter in the experiment, see e.g. Sect. 5.22 of [50]. 
Taking limiting values 

Q^i„ = 2GeV2, s,^in = 4GeV2, ymax = 0.85 (45) 

and 

|i|max = 0.5GeV2, |t,Uax = 0.9GeV2, (46) 

we find a Bethe-Heitler cross section of ubh = 1620 fb, which is between 15% and 20% smaller than 
the result we obtain in model 1 for ubh + crvcs + c^int with either an electron or positron beam. 
With an integrated luminosity of order lfb~ for HERMES running on a proton target [51], we 
deem this cross section to be too small, since it will be further decreased by experimental acceptance 
cuts and detection efficiency, and since according to (|35p the extraction of the beam charge or beam 
polarization asymmetry requires a differential measurement at least in the angle (f)-,^. Loosening the 
requirements (j45p or (j46p would increase the rate at the price of going to kinematics where the 
theoretical interpretation used in this paper becomes increasingly questionable. 

Higher luminosities than at HERMES can be achieved by the experiments at Jefferson Lab. With a 
currently available beam energy of up to E'g = 6 GeV, the requirements (145p and (j46p leave no available 
phase space, as can be seen from the bounds on x.,^ in (J14p and (J16p . This will be changed with the 
energy upgrade to -Eg = H GeV, which we consider in the remainder of this section. We assume that 
the outgoing electron, photon, and pion are detected experimentally. Note that identification of the 
pion (or of the recoiling neutron) is necessary to distinguish the signal process ep -^ e^jn^n from 
DVCS on the proton, ep -^ e'jp, which has a far greater rate. In the proton rest frame we have 






+ 2C0S(</)^ + Ve - V'n) ^(1 - X^)(l - X^){to - t)(t^o " t^) \ + o(^) 



cos^.= (l-i^^ii^il^) 



Q 

""^ +0(^), (47) 



VE^n-mi vq; 



where Et^ is the energy of the final-state pion and iDj^ the angle between its momentum and the lepton 
beam direction. Likewise, we find 

Q'^ , rnr \ d 1 2(l-2/)4"iIr -~^m- 
h Oim) , cosu^ = 1 -TT-^ — - 



E, = —^ + 0{m), cos^^ = l-^ n9 +^Vn3 

2xB ruN Q ^Q 

for the energy of the outgoing photon and its polar angle. To estimate the cross section and its 
dependence on the beam charge and polarization, we assume some minimal experimental cuts, which 
correspond to the acceptance planned for the CLAS-I— |- detector [52], 

S;>500MeV, 8°<T9e<45°, 

E^ > 100 MeV , 2° < t9^ < 40° , 

E^ > 200 MeV , 5° < i?^ < 135° . (49) 
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Figure 6: Scatter plots for the energies and polar angles of the final-state pion and photon in the 
proton rest frame. The plots are generated from the Bethe-Heitler cross section within the kinematics 
specified bv (|45]) . (f49l) . and pH) . 



Typical values of these quantities are shown in the scatter plots of Fig. [6l These plots have been 
generated using the Bethe-Heitler cross section, which dominates the total rate as we shall see shortly. 
We have imposed the kinematic requirements (|45p and 



0.3 GeV^ 



0.7 GeV^ 



(50) 



where compared with (146]) we have taken smaller |t|max and |t7r|max, so that the one-pion exchange and 
the leading-twist approximations are better fulfilled. We see that the pion has small to intermediate 
energy and covers a large angular region, whereas the photon is energetic and strongly focused in the 
beam direction. In the corresponding scatter plot for the outgoing electron kinematics (not shown 
here) we have 1.6 GeV < E'^ < 3.4 GeV and 13° < 'd^ < 25°, so that the experimental cuts on these 
quantities in (|49|) have no influence. 

Figure [7] shows corresponding scatter plots for other relevant variable pairs. In the first two panels 
we see the values of Q^, s-,^, and x^ at which the Compton process 7*7r — > ^-k can be probed with 
Eg, = 11 GeV. The last panel shows that typical values of the squared vrn invariant mass M^„ are 
between 1.3 GeV and 2.2 GeV , which unfortunately includes the region of nucleon resonances. We 
will come back to this point shortly. 

In addition to the contributions ubh, o"vcS) o'int of the different subprocesses and their interference 
to the integrated cross section, we evaluate the weighted differences 



(51) 



rtCOS (J}t^ 


= / (i(/>7r COS(/>7r 


'd(T{ei = +1) 

d4>TT 


da{e£ = - 


-1)1 


c-sin 0^ 


= / dcj)-,, sin(^^ 


'da{Pe = +1) 
d<P^ 


da {Pi = 


-1)1 



of cross sections for different beam charge or beam polarization. In the approximation given by (j29|) . 
(I32p . (I35p . they are respectively proportional to Re 7^^ and luiTCj^ in the interference term. In the 
Bjorken limit, the propagator factor P in ([33|) becomes (f)-,^ independent, so that S^^ and S^^ 
respectively give the coefficients of cos(/>7r and siiKJ).,^ in 27r dam^ / dcp^^ . With Eg = 11 GeV and the 
kinematics delineated by ([15]) and ([50]) . we find, however, a clear cp^^ dependence for P. For this 
reason, the interference term (|35p also contributes to the cross section integrated over (j).„, as seen in 
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Figure 7: Scatter plots for various kinematic quantities, generated in the same way as in Fig. [6l 
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Table 1: The contributions of the Bethe-Heitler and Compton processes and of their interference 
to the integrated cross section, as weh as the weighted cross sections defined in (fSTI) . Results are 
evaluated in model 1 for different kinematic constraints in addition to the cuts (j49p . The signs of 
(Tint and S^^ refer to an electron beam {ei = —1). Limiting values of Q^, s.,^, t, t^, and M^„ are 
given in units of GeV , and cross sections in units of fb. 



Qmin 


STrmin 


1 * 1 max 


1*77 1 max 


?/max 


M2 . 
Trnmm 


o"bh 


'^vcs 


O-JNT 


riCOS 4>Tr 




2 


4 


0.3 


0.7 


0.85 


— 


18.4 


0.88 


-0.18 


0.39 


7.57 


2 


4 


0.3 


0.7 


0.8 





5.12 


0.29 


-0.09 


0.17 


2.17 


2 


4 


0.3 


0.7 


0.9 





45.6 


1.86 


-0.27 


0.64 


17.9 


2 


4 


0.2 


0.7 


0.85 





0.41 


0.016 


-0.002 


0.004 


0.16 


2 


4 


0.5 


0.7 


0.85 





105 


6.52 


-2.32 


5.00 


46.2 


2.5 


4 


0.3 


0.7 


0.85 





2.55 


0.103 


-0.010 


0.018 


0.96 


2 


5 


0.3 


0.7 


0.85 





0.30 


0.013 


-0.003 


0.008 


0.12 


2 


4 


0.3 


0.5 


0.85 





16.2 


0.69 


-0.09 


0.18 


6.30 


2 


4 


0.3 


0.7 


0.85 


1.5 


13.4 


0.67 


-0.19 


0.42 


5.72 


2 


4 


0.3 


0.7 


0.85 


1.8 


5.08 


0.31 


-0.14 


0.30 


2.46 



Table [H We also recall that beyond the leading approximation in 1/Q, the weighted cross section 
5^™ receives a contribution from daycs in addition to the one from duiNT- This can be seen from 
the expressions in [25] and is well-known in the case of DVCS on a proton [22]. 

In Table [1] we show our results for the integrated and weighted cross sections calculated in model 
1 for different kinematic constraints. For the choice in (j45p and (j50p . shown in the first row, one 
obtains an integrated cross section of 18.4 fb. With an estimated luminosity of 3000 fb~ per year at 
CLAS-I--I-, this gives a very comfortable rate of about 55000 events, so that one may hope that even 
with realistic experimental cuts and detection efficiencies there will be sufficient statistics to perform 
a binning in several variables. 

The further entries in Table [T] illustrate how the situation changes if one modifies our baseline 
kinematic constraints. Raising ymax from 0.85 to 0.9 would more than double the rate, but as discussed 
above this would make the 1/Q expansion underlying the formulae ()32p and (|35p worse. Lowering 
instead ymax from 0.85 to 0.8, would reduce the rate by a factor of about 3.5. Increasing |t|max from 
0.3 GeV to 0.5 GeV we obtain a significantly higher rate, but the one-pion exchange approximation 
is much more problematic in that case. Decreasing |t|max to 0.2 GeV would improve the quality of 
the one-pion exchange approximation, but the resulting loss of rate is probably too much for such 
a cut to be useful. Taking a stricter cut of 2.5 GeV instead of 2 GeV for Q^j^ would make the 
leading-twist analysis of DVCS safer but decrease the cross section by about a factor of 7. An even 
stronger decrease is found if one requires s^r to be above 5 GeV instead of 4 GeV . By contrast, 
only very little rate is lost if one takes 0.5 GeV instead of 0.7 GeV for I^Trlmax; so that it may be 
worthwhile to consider a stronger cut on this variable. Finally, one may wish to impose a cut on 
the invariant mass of the vrn system, so as to reduce possible resonance effects in this channel, which 
are not taken into account in our theoretical description. The entries in the last two rows of Table [T] 
respectively correspond to a lower cut on M.„n at the mass of the A and at the mass of the A plus 
its width. The resulting rates show that at least part of the kinematic region where resonance effects 
may spoil a simple interpretation can be removed in an analysis without losing too much signal. 
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Table 2: As Table [T] but for different models of the pion GPDs. The kinematic constraints in ([^51) . 
(j49p . and ()50p are always assumed. The corresponding Bethe-Heitler cross sections is a^^ = 18.4 fb. 



model 


^vcs 


<7lNT 


qCOS (p^ 


r-sin 0^ 


1 


0.88 


-0.18 


0.39 


7.57 


2 


0.67 


-0.66 


1.57 


6.44 


3 


1.12 


-0.21 


0.43 


8.57 


4 


0.70 


0.25 


-0.62 


6.78 



In Fig. [8] we show the Bethe-Heitler cross section differential in t, t^r, or in x^. We see that in 
the kinematics of Jefferson Lab at 11 GeV a cutoff |t|max much below 0.3 GeV severely restricts the 
available phase space. The iyr spectrum falls off much more slowly than the corresponding t spectrum 
for DVCS on the proton, which readily follows from the slower decrease of the pion electromagnetic 
form factor compared with the proton form factors. Nevertheless, the falloff for |t7r|^0.2GeV is 
sufficiently steep to account for the weak dependence of the cross section on |t7r|inax seen in Table [TJ 
The spectrum in Xj^ is rather featureless and reflects the form of the pion flux factor in Fig. [2] together 
with the phase space boundaries Xmin and x^ax given in (fH|) and (fT6|) . 

The sensitivity to the pion GPDs of the Compton cross section and the interference term is 
documented in Table [2] for the different models introduced above. To obtain a more detailed picture, 
we plot in Fig. [9]the weighted cross sections (j5ip differential in t, t^r, and x^r. We see that the spread 
between models is much more pronounced for S™^ , with different signs and even a zero crossing in 
iyr for some of the models. By contrast, the variation of S^j^'^ is less drastic and concerns the size 
of the weighted cross section more than its dependence on t, i^r; or x^. On the other hand, the beam 
spin asymmetry S^j^'^ is significantly larger than the beam charge asymmetry S^^ and may be 
easier to measure in practice. 

6 Summary 

We have investigated the possibility to study DVCS on the pion in the reaction ep -^ e'yir^n. Such 
a study is experimentally demanding for several reasons. Firstly, the phase space is limited by the 
requirements of small t on one side (so that pion exchange dominates the process and the pion is not 
too far off-shell) and of large Q^ and Sj^ on the other side (so that an analysis based on the Bjorken 
limit is applicable). As can be seen from (J14p and (J16p . this favors experiments with a higher ep cm. 
energy. Secondly, the cross section for ep — > e^n^n is significantly smaller than the one for ep — > ejp, 
which puts high demands on both the luminosity and the experimental identification of the final state. 
We find that conditions for an experimental study of this process are not favorable in current 
experiments, with HERMES being limited by the available event rate and Jefferson Lab by the beam 
energy. After the planned energy upgrade to 11 GeV at Jefferson Lab, it should, however, be possible 
to investigate the reaction in detail. Using acceptance cuts relevant for the CLAS++ detector, we find 
comfortable event rates for kinematical conditions that may not be optimal but should be adequate for 
a first look at DVCS on the pion in the Bjorken regime and at intermediate values of the skewness ^. 
Using simple models for the GPDs of the pion, we estimate that information about them could be 
provided both by the beam spin and by the beam charge asymmetry, with the latter showing a more 
pronounced sensitivity to the ansatz for the GPDs but being smaller in size. Optimized studies might 
be feasible one day at a projected electron- ion collider [53], where in particular it should be possible 
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Figure 8: The Bethe-Heitler cross section taken differential in t, t^, or x^, and calculated with the 
cuts specified by (j45p . (09]), and (j50p . For the cross section differential in t we have omitted the first 
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Figure 9: As Fig. [8] but for the weighted cross sections (j5ip . calculated with the four GPD models 
specified in the text. The sign of S^^ refers to an electron beam. 
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to take the invariant mass of the vrn system above the resonance region. 

As a part of our model study, we have explored an ansatz for the pion GPDs that depends 
exponentially on t^r with a slope decreasing with x. Taking a functional form previously used for the 
nucleon [43] together with current parameterizations of the parton densities in the pion, we obtain 
an excellent fit to the experimental data of the electromagnetic pion form factor. Without adjusting 
further parameters, the ansatz then compares rather well with recent results from lattice QCD for 
the second Mellin moment of the GPDs. 
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